Search for: All records

Abstract Partition coefficients for rare earth elements (REEs) between apatite and basaltic melt were determined as a function of oxygen fugacity (fO2; iron-wüstite to hematite-magnetite buffers) at 1 bar and between 1110 and 1175 °C. Apatite-melt partitioning data for REE3+ (La, Sm, Gd, Lu) show near constant values at all experimental conditions, while bulk Eu becomes more incompatible (with an increasing negative anomaly) with decreasing fO2. Experiments define three apatite calibrations that can theoretically be used as redox sensors. The first, a XANES calibration that directly measures Eu valence in apatite, requires saturation at similar temperature-composition conditions to experiments and is defined by: ( E u 3 + ∑ E u ) Apatite  = 1 1 + 10 - 0.10 ± 0.01 × l o g ⁡ ( f o 2 ) - 1.63 ± 0.16 . The second technique involves analysis of Sm, Eu, and Gd in both apatite and coexisting basaltic melt (glass), and is defined by: ( Eu E u * ) D Sm × Gd = 1 1 + 10 - 0.15 ± 0.03 × l o g ⁡ ( f o 2 ) - 2.46 ± 0.41 . The third technique is based on the lattice strain model and also requires analysis of REE in both apatite and basalt. This calibration is defined by ( Eu E u * ) D lattice strain = 1 1 + 10 - 0.20 ± 0.03 × l o g ⁡ ( f o 2 ) - 3.03 ± 0.42 . The Eu valence-state partitioning techniques based on (Sm×Gd) and lattice strain are virtually indistinguishable, such that either methodology is valid. Application of any of these calibrations is best carried out in systems where both apatite and coexisting glass are present and in direct contact with one another. In holocrystalline rocks, whole rock analyses can be used as a guide to melt composition, but considerations and corrections must be made to either the lattice strain or Sm×Gd techniques to ensure that the effect of plagioclase crystallization either prior to or during apatite growth can be removed. Similarly, if the melt source has an inherited either a positive or negative Eu anomaly, appropriate corrections must also be made to lattice strain or Sm×Gd techniques that are based on whole rock analyses. This being the case, if apatite is primary and saturates from the parent melt early during the crystallization sequence, these corrections may be minimal. The partition coefficients for the REE between apatite and melt range from a maximum DEu3+ = 1.67 ± 0.25 (as determined by lattice strain) to DLu3+ = 0.69 ± 0.10. The REE partition coefficient pattern, as observed in the Onuma diagram, is in a fortuitous situation where the most compatible REE (Eu3+) is also the polyvalent element used to monitor fO2. These experiments provide a quantitative means of assessing Eu anomalies in apatite and how they be used to constrain the oxygen fugacity of silicate melts. 

The early Earth may have had oxidizing fluids feeding hydrothermal pools. 

Constraining the lithological diversity and tectonics of the earliest Earth is critical to understanding our planet’s evolution. Here we use detrital Jack Hills zircon (3.7 − 4.2 Ga) analyses coupled with new experimental partitioning data to model the silica content, Si+O isotopic composition, and trace element contents of their parent melts. Comparing our derived Jack Hills zircons’ parent melt Si+O isotopic compositions (−1.92 ≤ δ³⁰Si_NBS28 ≤ 0.53 ‰; 5.23 ≤ δ¹⁸O_VSMOW ≤ 9.00 ‰) to younger crustal lithologies, we conclude that the chemistry of the parent melts was influenced by the assimilation of terrigenous sediments, serpentinites, cherts, and silicified basalts, followed by igneous differentiation, leading to the formation of intermediate to felsic melts in the early Earth. Trace element measurements also show that the formational regime had an arc-like chemistry, implying the presence of mobile-lid tectonics in the Hadean. Finally, we propose that these continental-crust forming processes operated uniformly from 4.2 to at least 3.7 Ga.

 

The Vredefort impact structure, located in South Africa and formed 2.02 Ga, is the largest confirmed remnant impact crater on Earth. The widely accepted impactor diameter and velocity to form this crater are 15 km and 15 km/s, respectively, which produce a crater diameter of 172 km. This is much smaller than the most commonly cited estimates (250–280 km), and while previous results were able to match the geologic evidence known at that time, these impact parameters are not consistent with more recent geological constraints. Here, we conduct impact simulations to model the Vredefort crater formation with the shock physics code impact Simplified Arbitrary Lagrangian Eulerian (iSALE). Our numerical simulations show that combinations of the impactor diameter and impact velocity of 25 km and 15 km/s or 20 km and 25 km/s are able to recreate the larger crater size of ∼250 km. Moreover, these configurations can reproduce shock‐metamorphic features present in the impact structure today, including the distributions of breccia, shatter cones, planar deformation features in quartz and zircon, and melt. Our model also predicts that Vredefort and Karelia, Russia, where an ejecta layer from the impact was found, were approximately 2,000–2,500 km apart based on the layer thickness. Additionally, we use this model to predict the potential global effects of such a large impact by estimating the amount of climatically important gases released to the atmosphere at the time. Our work demonstrates the need to revisit previously estimated impactor parameters for large terrestrial craters in order to better characterize impact events on Earth and elsewhere.